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Abstract 


This  research  was  driven  by  the  need  of  the  Navy  to  improve  mine  countermeasures  in 
support  of  joint  littoral  warfare,  for  which  a  key  paradigm  is  to  locate  a  clear  or  low  mine-  and 
obstacle-density  path.  To  accomplish  this  goal,  it  is  imperative  that  one  be  able  to  predict  the 
behavior  of  mines  in  shallow  water  under  a  variety  of  coastal-ocean  conditions.  A  laboratory 
experimental  program,  coupled  with  theoretical  modeling,  was  conducted  at  Arizona  State 
University  for  the  purpose  of  improving  scientific  understanding  and  predictive  capabilities  of  the 
dynamics  of  the  motion  and  scour/burial  of  cobbles/mines  within  very  shallow  water  in  and  near 
the  surf  zone.  Results  of  these  studies  are  summarized  in  this  report. 

To  model  the  oscillatory  water  motion  in  different  regions  of  the  surf  zone,  three 
experimental  installations  were  constructed  and  used:  (i)  dam  break  tank  for  flow  in  a  swash 
region,  (ii)  standing  wave  tank  for  oscillatory  flow  near  the  surf  and  (iii)  large  wave  tank  to  model 
flow  in  the  entire  surf  zone,  including  wave  breaking  region.  The  data  on  the  motion  of 
cobbles/mines  along  solid  bottom  at  different  flow  conditions  were  collected  and  the  results  were 
explained  theoretically,  hi  the  case  of  a  sandy  bed,  detailed  data  on  the  bottom  morphodynamics 
(sand  ripples  and  bars  formation  and  evolution)  were  collected  and  explained  theoretically. 
Different  scenarios  of  cobble/mine  behavior  on  a  sandy  bed  were  simulated  experimentally  and 
explained  theoretically.  These  include:  (i)  slow,  mostly  onshore,  drift  of  cobbles,  (ii)  scour  around 
cobbles,  (iii)  periodic  burial  under  drifting  ripples,  and  (iv)  permanent  burial  under  migrating  sand 
bars. 
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1.  Introduction 

During  the  past  several  years,  the  Office  of  Naval  Research  (ONR)  has  supported  a 
research  program  at  Arizona  State  University  (ASU)  aimed  at  understanding  the  dynamics  of 
large  bottom  particles  (cobbles/mines)  subjected  to  periodic  background  flows  akin  to  those 
found  along  the  sea  floor  in  coastal  environments.  This  problem  has  important  applications  to 
mine  countermeasure  efforts  of  the  U.S.  Navy  in  littoral  zone.  Following  the  Gulf  War  (Desert 
Storm),  the  U.S.  naval  strategy  has  been  increasingly  focused  on  littoral  warfare,  with  major 
emphases  on  operational  and  naval  fleet  aids  in  coastal  environments.  Providing  tactical  support 
in  such  environments  requires  reliable  predictive  capabilities  of  the  fate  and  distribution  of 
mines,  which  can  be  accomplished  using  probabilistic  models  of  mine  behavior.  Development 
of  these  models  requires  a  fundamental  understanding  of  the  physics  and  dynamics  of  flow- 
mine-sand  interactions,  and  this  understanding  at  the  time  of  initiation  of  our  research  was 
insufficient. 

Taking  this  into  account,  a  comprehensive  research  program  was  directed  primarily  to 
understanding  the  (i)  dynamics  of  cobbles/mines  on  solid  bed  with  artificial  roughness,  (ii) 
morphodynamics  of  sandy  bed,  and  (iii)  cobbles/mines  scour/burial  mechanisms  on  a  sandy  bed. 
Parameterization  of  theses  processes  in  conditions  similar  to  those  that  occur  in  and  near  the  surf 
zone  was  also  of  interest.  The  focus  was  primarily  on  smaller  cobbles/mines  in  the  shallow  surf 
zone  (e.g.,  anti-tank  mines),  rather  than  large  “non-walking”  anti-ship  mines  in  deeper  coastal 
waters,  the  latter  being  the  theme  of  the  ONR  Mine  Burial  Program.  These  efforts  paralleled  the 
objectives  of  the  Organic  Mine  Countermeasures  (OMCM)  Future  Naval  Capabilities  (FNC) 
program  of  the  Office  of  Naval  Research. 

Over  the  years,  extensive  research  has  been  conducted  on  the  dynamics  of  wave  motion 
in  and  near  the  surf  zone  (see  the  schematic  in  Fig.  1)  and  on  the  resulting  sediment  transport 
(Louguet-Higgins  &  Stewart,  1964;  Yalin,  1977;  Engelund  &  Fredsoe,  1982;  Peregrina,  1983; 
Sleath,  1984;  Craik,  1985;  Battjes,  1988;  Mei,  1982,  1985,  Mei  &  Liu,  1993;  Mei  &  Yu,  1997; 
Fredsoe  &  Deigaard,  1992;  Nielsen,  1992;  Belorgey  et  al.,  1993;  Arcilla  et  al.,  1994;  Dean  & 
Dalrymple,  1994;  Yu  &  Mei,  2000;  Ridler  &  Sleath,  2000).  Most  of  these  studies  have 
concentrated  on  fine  sediment  transport,  but  only  little  attention  has  been  given  to  much  larger 
sediment  particles  of  the  size  of  cobbles.  (Cobbles  are  defined  by  the  American  Geophysical 
Union  as  sediment  particles  with  typical  diameter  D  =  6.4  -  25.6  cm;  Sleath,  1984).  Besides 
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obvious  engineering  and  geophysical  applications,  the  transport  of  such  large  particles  is  a  topic 
of  current  interest  to  the  Navy,  given  that  their  behavior  has  similarities  to  that  of  anti-tank  mines 
placed  in  coastal  waters  (Lott  &  Poeckert,  1996). 


Br*ik»r  zon« 


Shor* 

brtftk 


Reformed 

WIVM 


Swath 


Figure.  1.  Schematic  showing  different  regions  of  the  surf zone  ( Nielsen ,  1992). 


Motivated  by  the  important  naval  problem  of  being  able  to  predict  the  behavior  and  fate 
of  relatively  small  land  mines  in  and  near  the  surf  zone,  a  research  program  supported  by  the 
Coastal  Dynamics  Program  of  the  ONR  was  initiated  at  ASU.  In  contrast  to  the  rather  extensive 
literature  that  exists  on  the  transport  of  fine  sediments  in  the  coastal  region,  prior  to  the  said 
program  at  ASU,  very  little  research  had  been  done  on  the  behavior  of  cobbles/mines  on  sloping 
sandy  beaches.  During  the  period  of  this  grant,  the  ASU  group  made  steady  progress  in 
developing  a  better  understanding  of  cobble/mine  behavior  and  their  burial  as  driven  by  the  surf. 

The  long-ranee  goals  of  this  research  were  to  develop,  using  laboratory  experiments, 
theoretical  and,  in  some  cases,  numerical  analyses,  a  basic  understanding  and,  eventually,  a 
predictive  capability,  of  the  behavior  of  cobbles/mines  in  the  shoaling,  wave-breaking  and  swash 
regions.  Disk-shaped  anti-tank  mines  are  of  comparable  size  and  density. 

The  scientific  objectives  of  our  research  were  directed  toward  better  understanding  the 
behavior  of  cobbles  on  beaches,  which  are  permeable  with  movable  sand  bedforms.  These 
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objectives  permited  a  better  understanding  of  (i)  the  long-time  dynamics  of  bottom  topography 
(morphodynamics)  including  sand  ripples  and  bars  formed  in  a  model  surf  zone,  and  (ii)  the 
processes  related  to  scour  and  periodic/permanent  burial  of  cobbles  on  a  sandy  beach  with 
ripples  and  bars  in  time-  and  space-dependent  surf-zone  flow.  In  particular,  we  studied  (i)  the 
evolution  of  an  initially  flat  sandy  beach;  (ii)  the  long-time  behavior  of  the  bottom  topography, 
and  (iii)  the  behavior  of  model  cobbles  on  solid  and  sandy  beds 

The  approach  in  this  research  was  mostly  experimental.  Considering  the  nascent  nature 
of  research  on  this  problem,  a  step-by-step  approach  from  simplified  and  idealized  situations  to 
more  realistic,  and  subsequently  much  more  complicated  configurations,  was  adopted.  Realizing 
that  carefully  designed  experiments  are  an  avenue  for  achieving  the  scientific  goals,  a  number  of 
experimental  facilities  for  cobble/mine  research  were  designed  and  developed  at  Arizona  State 
University.  These  include:  (i)  three  different  experimental  facilities  capable  of  reproducing 
different  basic  background  flows  that  occur  in  and  near  the  surf  zone;  (ii)  a  system  of  devices  to 
control  and  measure  the  mean  and  turbulent  flow  characteristics  and  other  parameters;  and  (iii) 
instrumentation  for  employing  different  methods  of  flow  visualization  (see  Section  2  below). 

The  initial  studies  focused  on  cobble/mine  motion  along  impermeable  beaches.  At  first,  a 
series  of  experiments  was  conducted  in  which  the  dynamics  of  cobbles/mines  placed  on  a  solid 
slope  with  artificial  roughness  was  studied  in  background  flows  similar  to  those  observed  in 
different  regions  of  the  surf  zone.  Large  standing  waves  and  turbulent  bores  were  used  to  model, 
respectively,  steady  oscillations  and  swash  flows  (Luccio  et  al.,  1998;  Voropayev  et  al.,  1998). 
A  layer  of  sand  was  then  introduced  at  the  bottom  of  the  water  tank  and  the  dynamics  of  sand 
ripples  and  behavior  of  cobbles/mines  were  studied  under  oscillatory  flow  conditions 
(Voropayev  et  al.,  1999).  This  work  acted  as  a  precursor  to  the  studies  of  cobble/mine 
scour/burial  on  sloping  sandy  beaches. 

The  presence  of  sand  significantly  complicates  the  problem.  When  the  fluid  velocity 
exceeds  some  critical  value,  ripples  are  formed  on  an  initially  flat  sandy  bed,  which  affect  the 
behavior  and  burial/scouring  of  cobbles.  Although  the  physical  mechanism  of  ripple  formation 
is  qualitatively  well  understood  (see  Horikawa  &  Wanatabe,  1967,  Sleath,  1976;  1984,  1995  and 
references  therein),  owing  to  the  complexity  of  the  problem  a  quantitative  theory  in  this  regard 
was  still  in  infancy.  In  particular,  linear  and  weakly  nonlinear  analyses  have  been  conducted 
only  recently  (Blondeaux,  1990;  Vittory  &  Blondeaux,  1990;  Blondeaux  et  al.,  2000)  for  the  case 
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of  small  rolling-grain  ripples.  These  are  observed  at  small  velocities  of  fluid  oscillations.  At 
larger  velocities  the  height  of  the  rolling-grain  ripples  increases  and  flow  separation  occurs.  As  a 
result,  relatively  large  so-called  vortex  ripples  are  formed  (Fig.  2). 


Figure  2.  Large  vortices  shed 
periodically  from  symmetric  sand 
ripples  that  are  formed  in  symmetric 
oscillatory  flow  induced  by  standing 
waves.  The  water  moves  (a)  from  right 
to  left  and  (b)  from  left  to  right 
(Voropayev  et  al.  1999). 


Analytical  theories  for  such  complicated  nonlinear  flows  remain  untenable;  empirical  and 
semi-empirical  models  of  varying  complexities,  however,  are  in  frequent  use  (Sleath,  1984, 
1995;  Nielsen,  1992).  Needless  to  say,  the  presence  of  a  movable  sandy  bed  and  its  time 
evolution  extensively  modify  the  behavior  of  cobbles/mines.  Our  experiments  (Voropayev  et  al., 
1999)  show  that,  owing  to  large-time  instabilities,  the  bed  shape  does  not  reach  a  steady  state. 
The  ripples  migrate  with  time  and  their  migration  leads  to  periodic  burial/scouring  of  cobbles. 

More  realistic  modeling  of  typical  conditions  that  are  observed  in  the  oceanic  surf  zone 
required  large-scale  experiments.  Supported  by  ONR,  a  large  wave  tank  with  a  slope  was 
constructed  at  ASU.  A  series  of  experiments  was  conducted  and  theoretical  models  for 
cobble/mine  behavior,  which  were  developed  and  tested  previously  for  simplified  flow 
configurations,  were  tested  and  improved  for  surf  conditions  with  a  solid  bottom  (Voropayev  et 
al.,  2001).  Surf  conditions  with  a  sandy  slope  were  then  considered  (Voropayev  et  al.,  2000, 
2003).  Based  on  the  results  of  these  experiments,  a  physical  understanding  of  a  wide  range  of 
phenomena  was  attained  and  theoretical  models  describing  them  were  developed  and  tested.  The 
principal  results  of  these  studies  are  reported  in  our  publications  (which  are  given  at  the  end  of 
this  report)  and  summarized  in  the  next  sections. 

2.  Experimental  installations  and  flow  characteristics 
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Mostly  disc-shaped  cobbles/mines  were  used  in  the  experiments,  but  in  some  cases 
spherical  objects  were  also  used.  Because  the  spherical  cobbles  may  roll,  their  moment  of  inertia 
becomes  an  important  parameter  and  to  change  its  value  in  selected  experiments  spherical  shells 
filled  with  different  liquids  were  also  used.  In  total  more  than  20  different  cobbled  were  used  in 
experiments.  Typical  disc-shaped  cobbles  used  in  experiments  conducted  in  a  large  wave  tank 
are  shown  in  Fig.  3,  and  some  of  their  characteristics  are  given  in  Table  1  (diameter  D,  height  he, 
density  pc  and  static  Ks  and  dynamic  K  friction  coefficients).  Both  Ks  and  K  were  estimated  by 
dragging  with  relatively  small  velocity  the  submerged  cobbles  along  the  tank  floor,  and 
measuring  the  force  for  initiation  of  motion  and  force  in  a  steady  motion  by  a  standard  sensitive 
electronic  scale.  The  static  friction  coefficient  was  up  to  two  times  higher  than  the  dynamic 
friction  coefficient  (see  Table  1).  This  means  that  the  force  one  must  apply  to  initiate  the  motion 
is  twice  as  large  as  the  force  one  must  exert  to  keep  it  moving. 


Figure  3.  Model  cobbles/mines 
used  in  the  experiments.  Scale  is 
in  cm.  Other  parameters  are 
given  in  Table  1. 


Table  1 


D  (cm) 

h(cm) 

pc  (gem'3) 

K 

Ks 

1 

19.5 

8.5 

1.21 

0.5 

0.8 

2 

21.7 

9.1 

119 

0.3 

0.5 

A 

10.3 

5.3 

1.18 

0.5 

0.8 

B 

11.7 

4.7 

114 

0.5 

0.8 

C 

13.0 

5.4 

1.11 

0.5 

0.8 

To  model  flows  in  different  regions  of  a  surf  zone,  three  different  experimental 
installations  were  constructed  and  used  in  these  studies. 

Swash  flow  In  the  first  installation  a  dam  break  flow  was  used  to  model  flow  in  a  swash 
region  (Fig.  4).  In  this  approach,  a  dam  or  gate,  which,  prior  to  the  experiment,  blocks  a 
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reservoir  of  water,  is  nearly  instantaneously  removed.  This  results  in  a  well-defined  and 
reproducible  bore  which  moves  through  the  uniform  channel  of  rectangular  cross-section.  By 
varying  the  height  Ho  of  the  water  behind  the  dam,  the  height  and  propagation  speed  of  the  bore 
can  be  varied  (Fig.  5). 


(a)  (b) 

Figure  4.  Schematic  diagram  of  (a)  the  experimental  dam  break  flow  set-up 
and  (b)  model  cobble;  1-tank,  2-sloping  bottom,  3-cobble,  4-gate,  5-water. 


Two  similarities  can  be  noted  between  a  breaking  wave  impinging  on  a  beach  and  a 
turbulent  bore;  i.e.,  shape  and  speed.  Apart  from  these  large-scale  visual  similarities,  there  is 
also  an  internal  dynamical  similarity;  namely,  the  energy  lost  through  the  dissipation  in  a  bore 
gives  realistic  estimates  of  the  energy  lost  in  the  surf  zone  (e.g.,  see  Fredsoe  &  Deigaard,  1992). 
These  similarities  are  the  principal  reasons  the  dam-break  method  was  chosen  to  model  the  water 
motion  in  a  swash  zone. 


Figure  5.  Dam-break  flow:  Solid 
line  -  water  surface  without 
friction,  dashed  line  -  correction 
for  bottom  friction  at  the  tip  of  the 
front. 


The  great  advantage  of  this  method  is  that  the  background  velocity  produced  by  the  dam 
break  flow  is  reasonably  well  analyzed  theoretically  (Whitham,  1955,  1974)  and  to  a  good 
approximation  can  be  parameterized  as 

u  =  ai(gH0)1/2 
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(ai  =  const  ~  2/3,  see  details  in  Luccio  et  al.,  1998). 

In  the  experiments  the  cobbles  were  placed  at  a  fixed  distance  from  the  gate  along  the 
tank  axis  (Fig.  6a).  The  gate  was  opened  and  the  water  front  was  created.  The  distance  from  the 
gate  to  the  initial  position  of  the  cobble  was  the  same  (50  cm)  in  all  runs.  This  distance  was 
chosen  empirically  in  order  to  provide  a  tip  of  approximately  constant  height  as  the  water  front 
reached  the  cobble.  Upon  impingement  of  the  water  front,  (Fig.  6b),  the  cobble  started  to  move 
(Fig.  6c)  with  a  velocity  U  =  dxc/dt,  where  xc  is  the  cobble  displacement  from  the  initial 

position  at  t  =  0. 


Figure  6.  The  motion  of  a  cobble  in  a 
dam-break  flow:  (a)  -  top  view,  (b-c)  - 
side  view;  (a)  - 1  <  0,  (b)  - 1  =  0,  (c)  - 1  = 
0.04  s,  (d)  -  t  =  0.67  s.  The  flow  moves 
from  right  to  left.  The  time  origin,  t  =  0, 
is  chosen  at  the  moment  when  cobble, 
which  was  initially  at  rest  at  x  -  50  cm, 
starts  to  move.  Experimental  conditions: 
H0  =  20  cm,  D  =  10  cm,  hc  =  4  cm,  pc  = 

1.5  gm  cm~3. 


Two  separate  tank  floors  were  used  in  the  experiments;  i.e.,  smooth  (Plexiglas)  and 
rough,  for  which  a  thin  layer  of  sand  was  glued  to  the  Plexiglas.  This  permitted  a  variation  of  the 
kinetic  friction  coefficient  K  between  the  bottom  and  cobble  by  a  factor  of  four.  By  tilting  the 
bottom  at  a  small  angle  to  the  horizontal,  the  effect  of  a  sloping  bottom  was  modeled  in  the 
experiments.  A  total  more  than  40  experiments  were  conducted  with  different  initial  conditions 
and  data  on  the  cobble  displacements  as  a  function  of  time  were  taken  using  a  video  camera  with 
high  resolution. 

Oscillatory  flow  was  modeled  using  standing  waves  induced  in  a  large  (length  L  =  366 
cm)  tank  of  rectangular  cross-section  (Fig.  7).  The  waves  were  forced  by  a  vertical  paddle 
installed  at  one  end  of  the  tank  and  connected  through  a  mechanical  drive  to  a  precision  motor; 
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the  paddle  oscillates  with  a  prescribed  amplitude  and  frequency,  CD,  thus  generating  periodic 
progressive  surface  waves  of  the  amplitude  a/2.  When  these  waves  reach  the  end  wall  of  the 
tank,  they  reflect  without  a  significant  loss  of  amplitude  and  propagate  in  the  opposite  direction 
with  the  same  frequency.  The  superposition  of  these  progressive  waves  leads  to  a  standing  wave 
pattern.  The  vertical  velocity  near  the  central  portion  of  the  tank  is  small,  while  the  horizontal 
velocity  u  is  practically  depth  independent  and  can  be  approximated  as  (Whitham,  1974) 

.  7CX  .  . 

u  «  ecosin — sin  cot «  ecosin  cot , 

L 

where  e  =  aL/(7th0)  (ho  is  the  water  depth).  Employing  standing  waves  of  the  first  mode,  it  was 
possible  to  produce  horizontal  oscillations  in  the  central  potion  of  the  tank  with  velocity 
amplitude  up  to  80  cms  \ 


Figure.  7.  Standing  wave  tank:  1 
-  tank,  2  -  water,  3  -  paddle,  4  - 
eccentric  drive  mechanism,  5  - 
solid  bottom  with  artificial 
roughness  or  layer  of  sand,  6  - 
cobble,  7  -  ADV probe. 


In  the  experiments,  a  cobble  was  placed  on  the  floor  in  the  vicinity  of  the  tank  center  and 
its  subsequent  motion  was  investigated.  Two  different  cobble  shapes  were  considered;  i.e, 
spheres  of  diameter  D  and  density  pc  and  disks  of  diameter  D,  height  he  and  density  pc.  Typical 

-3 

parameter  values  used  were  D  =  1.2-22  cm,  hc=  2-5  cm,  pc  =  1.04-7.9  gem  (for  details  see 
Voropayev  et  al.,  1998).  The  central  portion  of  the  tank  floor  was  removable  and  was  either 
smooth  or  rough  (sand-covered  plate).  This  allowed  a  variation  of  the  coefficient  of  friction 
between  the  cobble  and  the  floor  in  a  factor  of  5. 

The  model  cobble  (sphere  or  disk)  was  placed  on  the  solid  floor  near  the  center  of  the 
tank  and  its  position  (and  the  angular  velocity  for  spherical  cobbles)  was  monitored.  Then  a 
layer  of  sand  (depth  15  cm)  was  placed  at  the  bottom.  In  these  experiments  the  time  evolution  of 
the  sandy  bed  (without  cobbles)  was  studied  and  data  on  the  formation  and  drift  of  vortex  ripples 
were  collected.  Then  selected  cobbles  were  placed  at  sandy  bed  and  their  behavior  (motion, 


< - - -  L  - - - ► 
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scour/burial)  was  studied.  In  total  more  than  80  experiments  were  conducted  using  this  standing 
wave  tank 

A  large  wave  tank  was  used  to  model  the  water  motion  in  the  entire  surf  zone,  which 
includes  shoaling,  wave-breaking  and  swash  regions.  This  glass-walled  wave  tank  (length  32  m, 
width  0.9  m,  depth  1.8  m)  was  designed  and  constructed  at  ASU.  An  oblique  view  of  the  tank  is 
given  in  Fig.  8  and  the  schematic  of  the  experimental  system  is  shown  in  Fig.  9.  The  plane 
vertical  paddle,  installed  at  one  end  of  the  tank  and  driving  by  a  piston,  is  used  to  generate 
periodic  waves.  The  wave  maker  is  driven  by  a  hydraulic  system,  which  is  capable  of  producing 
a  periodic  sinusoidal  or  step-like  motion  with  a  highly  stable  frequency.  Periodic  sinusoidal 
waves  of  frequency  to  =  0.2-0.7  Hz  and  peak-to-peak  amplitude  of  the  paddle  motion  up  to  2e  — 
25  cm  were  used  in  most  of  the  experiments.  A  wooden  bottom  with  slope  1:24  (length'  26  m)  is 
installed  opposite  to  the  paddle.  In  experiments  with  solid  bottom  a  thin  layer  of  sand  is  glued 
on  top  of  the  sloping  bottom  to  increase  the  bottom  friction.  In  experiments  with  sandy  bed  a 
layer  of  sand  (thickness  23-25  cm)  was  placed  along  the  bottom.  . . 


Figure  8.  The  Wave  tank  facility  that  was  used  to  model  the  flow  in  the  surf  zone.  Numbers 
show  the  sections  as  counted  from  the  wave  maker,  which  is  at  the  right  end  of  the  tank.  The 
distance  between  vertical  bars  is  61  cm. 

A  typical  nonlinear  wave  propagating  from  right  to  left  along  sloping  bottom  is  shown  in 
Fig.  10.  In  contrast  to  dam  break  flow  (Figs.  4,5)  and  oscillatory  flow  induced  by  large  standing 
waves  (Fig.  7),  where  the  background  flow  velocity  can  be  estimated  theoretically,  the  flow 
characteristics  in  large  wave  tank  cannot  be  estimated  accurately  theoretically. 


i? 
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Figure  9.  Progressive  wave  tank  experimental  setup:  1  -  wave  tank,  2  -  water,  3  - 
vertical  wave-maker  (frequency  (0,  amplitude  of  horizontal  displacement  £q),  4  - 
slopping  solid  or  sand  bottom,  5  -  cobble,  6  -  Acoustic  Doppler  Velocimeter  (ADV) 
and  wave  gauge  attached  to  carriage,  7  -  hydraulic  system  to  move  the  wave  maker. 
Section  numbers  are  also  shown  and  the  length  of  one  section  is  61  cm. 


Figure  10.  A  typical  nonlinear  wave 
propagating  from  right  to  left  along 
a  slope  and  then  breaking  near  the 
break  point.  Two  model  mines  are 
seen  at  the  bottom  in  section  #29. 


Significant  research  has  been  reported  on  weakly  nonlinear  waves  in  wave  tanks 
(Fontanet,  1961;  Madsen,  1971;  Flick  &  Guza,  1980;  Svendsen,  1985),  but  until  now  there  is  no 
satisfactory  theory  to  accurately  predict  the  wave  characteristics  and  water  velocities  associated 
with  strongly  nonlinear  shoaling  waves  on  a  slope.  Although  the  wave-maker  forcing  is 
sinusoidal,  the  waves  generated  become  increasingly  nonlinear  as  they  propagate  along  the  slope 
and  increase  in  wave  steepness.  In  a  model  surf  zone  (Fig.  9)  the  background  flow  is  a  rather 
complicated  function  of  time  and  space  (see  Figs.  1 1,  12  below),  and  the  water  particle  velocities 
u(x,t)  associated  with  the  motion  of  strongly  nonlinear  waves  (Fig.  10)  propagating  and  breaking 
along  the  slope  is  rather  complicated.  Only  recently  some  progress  in  direct  numerical 
simulations  was  achieved  for  such  flows  and  u(x,t)  was  accurately  calculated  but  this  is  a  rather 
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complicated  procedure,  which  requires  accurate  preliminary  calibrations  (Grilli  at  al.,  2003).  As 
accurate  data  on  u(x,t)  are  needed  to  estimate  the  background  flow  characteristics,  direct 
measurements  were  made.  The  velocity  components  were  measured  using  a  standard  three- 
component  Acoustic  Doppler  Velocimeter  (ADV).  This  instrument  had  a  relatively  low  noise 

level  O  0.2  cms'1)  and  a  small  cylindrical  (diameter  0.6  cm,  height  0.3  cm)  measurement 
volume  (=  0.1  cm3)  with  a  data  collection  frequency  up  to  25  Hz  (for  detail  see  Snyder  & 
Castro,  1998).  The  wave  elevations  were  measured  using  standard  wave  gauges.  The  ADV 
probe  is  attached  to  a  long  supporting  rod,  which  is  fixed  to  a  carriage.  The  carriage  is  mounted 
on  rails  on  top  of  the  tank  and  has  two  degrees  of  freedom,  traversing  in  the  horizontal  and 
vertical  directions.  In  order  to  measure  the  Lagrangian  displacement  field  in  the  water,  small 
neutrally  buoyant  particles  were  introduced  into  the  flow.  To  obtain  representative  data,  the 
velocity  and  wave  elevation  records  were  phase  averaged  over  50-60  periods  of  oscillations.  An 
example  of  such  velocity  measurements  is  shown  in  Fig.  11. 


Figure  11.  Typical  horizontal,  u,  and 
vertical,  w,  components  of  the  water 
velocity  under  a  periodic  nonlinear  wave 
along  the  slope.  Section  #25,  depth  is  10 
cm  above  the  bottom.  The  data  are 
averaged  over  50  wave  periods  and  the 
wave  period  is  2.5  s. 


Horizontal  velocity  component  u  in  the  along  flow  direction  x  and  the  vertical  component 
w  in  the  vertical  direction  z  are  shown  in  Fig.  11.  There  is  practically  no  water  motion  in  the 
across  flow  direction  y  and  this  velocity  component  is  not  shown.  As  can  be  seen,  although  the 
paddle  motion  is  purely  sinusoidal,  the  water  particle  velocity  components  u  and  w  become 
strongly  asymmetric  as  the  wave  propagates  along  the  slope.  The  measurements,  similar  to  those 
shown  in  Fig.  1 1  were  made  in  different  sections  of  the  tank  and  typical  results  for  the  horizontal 
along-slope  water  velocities  u  in  sections  19-27  are  shown  in  Fig.  12.  As  can  be  seen,  the 
amplitude  of  the  horizontal  water  velocity  changes  along  the  slope.  The  maximum  velocity 
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increases,  while  the  time  interval  in  which  the  velocity  is  positive  during  each  period  decreases. 
This  is  the  result  of  the  wave  steepening  and  related  nonlinear  effects. 


As  can  be  seen  from  the  results  of  measurements  presented  above,  the  water  motion  in 
the  tank  is  strongly  asymmetric  (in  the  onshore  and  offshore  directions)  and  varies  with  the 
distance  along  the  slope.  These  are  the  most  important  differences  between  the  present  case  of 
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Figure  12.  Typical  horizontal  water 
particle  velocities,  u,  as  function  of  position 
in  the  wave  tank.  The  data  were  taken  in  |  20 

different  sections  at  depth  10  cm  from  the  =  o 

bottom.  Standard  phase  averaging  was 
applied.  The  waves  break  in  section  #30.  20 

Section  numbers  are  shown  in  the  legend.  ^ 

nonlinear  progressive  shallow  water  waves  propagating  along  a  slope  and  much  more  simple 

cases  of  swash  flow  generated  by  a  dam  break  method  and  of  spatially  independent  oscillatory 

flow  generated  by  large  standing  waves. 

The  data,  similar  to  those  shown  in  Figs.  1 1, 12,  were  digitized  and  used  as  inputs  for  the 
model  calculations,  which  are  discussed  below. 


3.  Cobbles/mines  on  impermeable  bed 

The  results  of  experiments  on  the  motion  of  cobbles/mines  along  horizontal  and  slopping 
solid  beds  for  different  background  flows  are  now  discussed.  The  general  problem  here  may  be 
formulated  as  follows.  Let  a  cobble/mine  rests  on  a  solid  bed  and  be  subject  to  the  action  of  (in 
general  arbitrary)  background  flow.  When  the  flow  velocity  exceeds  some  critical  value  the  object 
starts  moving.  What  are  the  critical  conditions  of  the  initiation  of  motion  and  how  the  object  will 
move?  The  partial  answers  to  these  questions  are  given  below. 

Based  on  the  results  of  numerous  experiments  conducted  under  simplified  background 
flow  conditions  (swash  flow  and  symmetric  oscillatory  flow)  theoretical  models  were  developed 
to  describe  the  motion  of  heavy  disks,  spheres,  spherical  shells  filled  with  liquid  (to  change  the 
moment  of  inertia)  in  these  basic  time-dependent  flows  (Luccio  et  al.,  1998;  Voropayev  et  al., 
1998).  Using  the  proper  theoretical  parameterizations  for  the  velocity  in  such  flows,  analytical 
solutions  were  derived  in  closed  form  and  verified  experimentally. 
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The  models  developed  for  these  configurations  are  rather  general  and  permit  one  to 
consider  more  complicated  cases  for  which  the  background  flow  depends  not  only  on  time  but  on 
the  spatial  (along-slope)  coordinate  as  well;  i.e.,  the  situation  in  a  real  surf  zone  or  that  modeled 
in  the  large  wave  tank  (Figs.9,  10).  In  cases  for  which  the  background  flow  is  primarily  a 
function  of  time  only,  it  was  possible  to  derive  analytical  solutions  describing  the  cobble/mine 
displacement.  For  the  general  case  in  which  the  background  flow  is  a  function  of  both  time  and 
the  spatial  coordinate  (Figs.  11,12),  analytical  solutions  were  untenable  and  the  resulting 
nonlinear,  integro-differential  equations  for  the  cobble/mine  motion  were  solved  numerically. 

Consider  briefly  a  general  case  in  which  the  background  flow  is  a  function  of  both  time 
and  the  spatial  coordinate  (for  details  see  Voropayev  et  al.,  2001).  A  schematic  of  the  system 
under  consideration  is  shown  in  Fig.  13.  Here  u(x,t)  is  the  background  water  velocity  in  the 
stream-wise  or  x-direction,  U(t)  cobble  velocity.  Fa  drag  force,  Fl  lift  force,  Ff  dynamic  bottom 
friction,  Fp  accelerating  pressure  term,  G  =  (Mc-M)g  reduced  gravity  force,  g  gravitational 
acceleration,  Me  mass  of  cobble,  M  water  mass  displaced  by  cobble,  a  slope  angle,  Fg  =  G  sina 
and  Fn  =  G  cosa. 


Figure  13.  Definition  sketch  with 
coordinate  system  and  periodic 
forces  that  act  on  a  cobble/mine 
placed  on  a  sloping  beach  in  a 
surf  zone. 


K-D-W 


Consider  a  disk-shaped  cobble  (Fig.  13)  of  diameter  D,  thickness  h  and  density  pc 
initially  (t  <  0)  at  rest  at  xc  =  0,  on  a  plane  sloping  at  an  angle  a  to  the  horizontal.  At  t  =  0,  the 
surrounding  fluid  of  density  p  (<  pc)  begins  to  move  with  an  arbitrary  velocity  u(x,t)  parallel  to 
the  floor.  If  the  net  force  that  acts  on  the  cobble  is  not  equal  to  zero,  the  cobble  begins  to  move 
with  a  velocity  U  =  dxc/dt  in  the  direction  of  the  x-axis. 

The  balance  of  momentum  for  a  cobble  in  the  stream-wise  or  x-direction  gives 

(Mc+KmM)^=(l  +  Km)M^-  +  Fd-Gsina  +  Ff,  (1) 

dt  dt 

(Km  is  the  virtual  mass  coefficient,  du/dt  =  du/dt  +  u3u/3x).  Using  standard  parameterizations 
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1  ,  ,  u  f.l,  Gcosa>Er,l 

Fd=>cd|u-U|(u-U),  F^-KEfGcosa-^,  B-|ft  Gcosa£^}. 


Ji=IpS2CL(u-U)2 


JiD2h 


7tD2h 


Mc=Pc— >M  =  p— 


(Si  =  Dh,  S2=7tD2/4,  Cd  and  Cl  are  the  drag  and  lift  coefficient,  K  is  the  kinetic  friction 
coefficient)  and  denoting  p*  =  pc/p,  A  =  u  -  U,  which  is  the  velocity  of  the  cobble  relative  to  the 
water  velocity,  and  using  the  identity  dA/dt  =  dA/dt  +  UdA/ dx ,  (1)  takes  the  form 

(p*  +  Km)(^ + Aji)  +  ^  A  |A|  =  (p*  - 1)  ^- + (p*  -  l)g  sin  a 

dt  dx  tcD  1  1  dt  ^ 

+  KB|fj<p'-1)gc°sa-^A2} 

For  the  configuration  considered,  Cl  ~  C^  «  Co  *  0.6  and  Km  =1.2  can  be  used  (Yalin,  1977; 

Sarpkaya,  1986;  Voropayev  et  al.,  1998, 2001). 

In  non-dimensional  variables 

*  u  tt*  U  A  *  ,  . 

u  = — ,  U  = — ,  A  = — ,  t  =  cot, 

eco  eco  ecu  ^ 

D*=— ,  h*=i  P*=££,  g‘=-^, 

e  e  p  eco2 

(the  amplitude  £  and  period  1/co  of  the  paddle  oscillations  are  used  as  the  length  and  time  scales) 
one  obtains 

.  *  T,  wdA  A*du*.  .  *  du  2C0  *  .*|  .  *  - v  ♦  . 

(P  +KJ(— +A— )=(p  -!)— r — rr  A  A  +(p  -l)g  sin  a 
dt  dx  *  dt  JtD  1 

r  ci  ’  w 

+  KER|(p*  -l)g*cosa - r(A*)2| 


where 


1,  (p*  -  l)g*  cos  a  >  -%  A*2 

,  2h 

0,  (p*-l)g*cosa<-%A*2 

I  2h 


U*  u*  -  A* 


When  u  does  not  depend  on  the  spatial  coordinate  (du/dx  =  0),  the  solution  for  A*  may 


be  found  analytically  and  the  non-dimensional  cobble  displacement  xc  =  xc  /£  is  given  by 
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(6) 


t*  t*  t* 

xc  =  JU  dt  =/u  dt  -  jAdt  . 

ooo 

Some  examples  of  practical  interest  are  given  in  Luccio  et  al.  (1998)  and  Voropayev  et  al. 
(1998). 

If  the  form  of  background  velocity  u(x,t)  is  known,  for  example,  from  direct 
measurements  as  in  our  case,  then  the  system  of  equations  may  be  closed.  Choosing  the  initial 
conditions,  for  example,  as 

A*  =  A*  (0)  and  x*  =  x*c  =  x*c(0)  at  t*  =  0, 

the  governing  equation  (4)  can  be  solved  using  a  fourth  order  standard  Runge-Kutta  method  and 
cobble  positions  may  be  calculated.  For  this  purpose  digitized  experimental  data,  similar  to 
those  shown  in  Figs.  11, 12,  were  used  as  inputs  for  u(x,t)  in  the  model  calculations. 

In  the  experiments  cobbles/mines  of  different  sizes  were  placed  along  the  floor  in  a  wave 
tank  (Fig.  8)  with  model  surf  flow  and  their  evolution  with  time  was  studied  and  compared  with 
the  model  predictions.  Onshore  and  offshore  mean  motions  of  cobbles  as  well  as  steady 
oscillations  with  zero  mean  displacement  were  observed  for  different  conditions  and  two  typical 
examples  are  shown  in  Figs.  14, 15. 


Figure  14.  Averaged  over  one  wave 
period  cobble  position  xc  as  a 
function  of  time  t  for  a  medium-size 
cobble  (cobble  B  in  Table  1). 
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Figure  15.  Mine  position  xc  as  a 
Junction  of  time  t  for  a  large 
model  mine  (mine  2  in  Table  1) 
provided  by  NRL.  Experimental 
data  were  taken  using  high- 
resolution  digital  imaging. 


«(•> 

For  the  range  of  parameters  used  in  the  experiments,  satisfactory  agreement  between  the 
measured  and  calculated  values  of  the  cobble/mine  displacements  as  a  function  of  time  was 
obtained  in  other  experiments  including  case  when  one  of  the  cobbles  was  lifted  (at  least  partly 
from  the  bottom  as  it  moves  forward  in  section  #29. 


Figure  16.  A  photograph  showing  the 
separation  of  cobble  2  ( replica  of  anti¬ 
tank  mine)  from  bed  during  onshore 
(from  right  to  left)  stage  of  motion. 


Overall,  the  behavior  of  cobbles  on  a  slope  depends  on  the  cobble  size,  density  and  initial 
cobble  position.  For  the  wave  parameters  used  in  the  experiments,  the  cobble/mine  does  not 
move  at  all  when  its  density  is  relatively  large  (pc  >  1.5  gem'3).  In  this  case  the  (static)  friction 
force  exceeds  the  net  driving  force  from  the  surrounding  fluid  at  any  position  of  cobble  on  the 
slope.  We  could  not  increase  significantly  the  driving  force  from  the  surrounding  fluid  (wave 
amplitudes)  without  the  corresponding  loses  of  wave  stability  and  for  this  reason  cobbles  with 
smaller  densities  (see  Table  1)  were  used  in  these  experiments.  These  lighter  cobbles  can  move 
along  the  slope,  mostly  onshore,  but  in  general  their  direction  of  motion  may  depend  on  the 
initial  cobble  position  and  cobble  parameters  (see  details  in  Voropayev  et  al.,  2001). 

4.  Sandy  bed  morphodynamics 

In  these  studies  a  layer  of  sand  (thickness  20-25  cm)  was  placed  on  the  bed.  Two  kinds 
of  quartz  sand  (density  ps  =2.6  gem'3,  mean  grain  size  d  =  0.04  cm  and  0.1  cm)  were  used  in  the 
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experiments.  The  sand  surface  was  leveled  and  the  background  flow  was  initiated.  Particular 
attention  was  paid  to  the  instability  of  the  sand-water  interface,  formation  of  sand  ripples  and  bar 
and  temporal  evolution  of  the  bottom  morphology  at  different  background  flow  conditions.  Two 
series  of  experiments  were  done.  The  first  series  was  conducted  in  symmetric  oscillatory  flow 
using  standing  wave  tank  (Fig.  7)  with  horizontal  bottom.  Second  series  of  experiments  was 
done  in  a  model  surf  flow  reproduced  in  a  large  wave  tank  with  a  slope  (Fig.  8). 

The  general  problem  here  may  be  formulated  as  follows.  Let  initially  flat  layer  of  sand  (on 
a  horizontal  or  sloping  bottom)  is  subjected  to  the  action  of  a  periodic  background  flow.  It  is  well 
known  that  at  some  critical  conditions  large  vortex  ripples  and  bars  are  formed  at  the  sand  surface. 
How  fast  these  ripples  and  bars  are  formed,  what  are  their  quasi-steady  characteristics,  their 
migration  velocity  in  different  flow  conditions?  Partial  answers  on  these  questions  are  given 
below/ 

Qualitative  observations  of  the  experiments  can  be  summarized  as  follows.  After  the 
waves  are  generated,  an  oscillatory  shear  stress  is  developed  near  initially  flat  bottom  (Fig.  17a) 
and  rolling  grain  sediment  transport  was  observed  when  the  velocity  amplitude  exceeds  a  critical 
value.  Owing  to  the  instability  at  the  water-sand  interface,  small  (typical  height  0.5  cm)  three- 
dimensional  ripples  (typical  horizontal  spacing  of  4-5  cm)  are  generated  at  the  bed.  With  time, 
the  height  of  the  initial  3D  ripples  increases;  the  ripples  then  merge,  forming  the  well-known  2D 
vortex  ripples  oriented  across  the  flow  (Fig.  17b). 


Figure  17.  Side  view  photographs 
showing  (a)  •  initially  flat  sandy 
bottom;  (b)  -  growing  and  (c)- 
quasi- steady  sand  ripples  formed 
on  a  slope  under  progressive 
waves  with  relatively  large 
intensity..  The  horizontal  size  of 
frames,  61  cm,  gives  a  length 
scale. 
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These  2D  ripples  initially  have  a  relatively  small  height  h(t)  and  spacing  L(t),  but  with 
time  their  size  increases  and,  after  an  adjustment  time  interval  (see  below),  they  reach  a  quasi¬ 
steady  equilibrium  state  (Fig.  17c)  with  a  characteristic  height  h0  and  spacing  Lo.  With 
increasing  ripple  size  and  height  to  length  ratio,  h/L,  the  flow  around  the  ripples  becomes 
turbulent  and  intense,  and  vortices  are  periodically  generated  near  the  ripple  crests.  Ripples  in 
this  stage  are  called  the  “vortex  ripples.”  When  the  waves  are  symmetric  (relative  to  the  vertical 
plane),  as  in  the  case  of  large  symmetric  standing  waves,  vortices  are  formed  periodically  and 
symmetrically  at  both  sides  of  the  ripples  which  also  tend  to  be  symmetric  (Fig.  2).  Similar  (but 
slightly  asymmetric)  ripples  are  also  observed  for  progressive  waves  on  a  slope  far  offshore  from 
the  break  point  where  the  wave  asymmetry  is  relatively  small.  When  the  wave  asymmetry 
increases,  the  vortices  are  formed  mostly  at  the  onshore  side  of  the  ripple  while  they  do  not  occur 
at  the  offshore  side;  the  resulting  ripples  are  strongly  asymmetric  (Fig.  17c). 

Note  that  for  symmetric  flow  ripple  generation  occurs  uniformly  along  the  horizontal 
bed.  For  progressive  waves  on  a  slope,  at  early  times  3D  ripples  are  generated  rapidly  just 
offshore  of  the  breaking  point  were  the  velocity  amplitude  has  its  maximum.  The  ripple- 
formation  zone  then  propagates  as  a  front  in  the  offshore  (down-slope)  direction;  after  some  time 
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the  entire  slope  is  covered  with  ripples  (Fig.  18).  ; 


Figure  18.  Ripple  front  (shown 
by  arrow)  propagating  (b,  c) 
along  initially  flat  slope  (a)  in 
the  offshore  (from  left  to  right) 
direction  under  progressive 
waves  on  a  slop. 


The  front  propagation  velocity  strongly  depends  on  the  wave  characteristics  and  water 
depth  and  decreases  with  the  distance  X  from  the  break  point.  If  the  waves  are  not  very 
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energetic,  than  after  some  typical  time,  to,  the  ripple  front  stops  propagating  along  the  slope  at  a 
distance  So  from  the  point  of  ripple  formation  near  the  break  point.  As  a  result,  part  of  the  slope 
is  covered  with  ripples  while  the  rest  remains  flat  for  a  long  time  (4-5  hours).  For  very  energetic 
waves,  Xo  is  rather  large  and  ripple  front  propagates  along  the  entire  slope.  Experimental  data  on 
the  front  propagation  were  collected  and  explained  using  a  model  developed  by  the  present 
investigators.  The  underlying  concept  of  the  model  is  that  the  ripple  height  h  increases  with  time 
in  accordance  with  (7)  (see  below)  while  characteristic  time  x  vary  along  the  slope  in  accordance 
with  the  background  flow  variation  (ho  increases  and  x  decreases  in  the  offshore  direction).  As  a 
result,  the  non-dimensional  front  distance  increases  with  time  as 

X/Xo  =  1  -  exp(-t/x), 

where  x  is  a  characteristic  time  that  depends  on  the  wave  frequency  and  mobility  parameter  (see 
details  in  Voropayev  et  al.,  2003). 

Analyses  of  experimental  data  (see  Fig.  19)  show  that  in  both  cases  (horizontal  bed  and 
sloping  bed)  the  non-dimensional  length  L/Lo  and  height  h/ho  of  the  generated  ripples  increase 
with  the  non-dimensional  time  t lx  as 

L/L0  =  1  -  exp(-  t/x) ,  h/h0  =  1  -  exp(-  t/x) ,  (7) 

where  L  and  h  are  the  ripple  current  length  and  height,  and  Lo  and  h0  are  their  equilibrium  quasi¬ 
steady  values,  which  may  be  calculated  using  well  known  empirical  formulas  (see,  e.g.,  Sleath, 
1984, 0'Donoghue  &  Clubb,  2001)  and 

X  =  C/cD\)/1/2 

is  a  characteristic  time  scale  for  ripple  formation  introduced  in  Voropayev  et  al.  (1999,  2003). 
Here  ©  is  the  wave  frequency,  C  =  2500  is  an  empirical  constant  and  \| f  is  the  mobility 
parameter,  defined  as 

V  =  U2/g(s-l)d  (8) 

(s  =  ps/p ,  d  is  the  mean  sediment  diameter,  g  is  the  gravitational  acceleration,  U  is  the  amplitude 
of  fluid  velocity). 

The  main  difference  between  two  cases  is  that  for  a  horizontal  bed  (Fig  19a)  \\f  doesn’t 
depend  on  the  position  x  along  the  bed,  while  for  a  slope  (Fig.  19b)  \|/  increases  along  the  slope 
in  the  onshore  direction  (because  the  maximum  along  flow  velocity  U  increases). 
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Figure  19.  The  non-dimensional  ripple  lengths  L/Lo  as  functions  of  the  non-dimensional  time  t/rfor 
different  experiments  (symbols)  conducted  with  (a)  oscillatory  flow  on  horizontal  bed,  (b)with 
progressive  waves  along  the  slope.  The  solid  line  shows  the  approximation  function  (7).  In  (b)  the 
data  were  taken  at  different  incoming  wave  characteristics  and  in  different  sections  along  the  slope  and 
changes  in  the  values  of  ynv  (8)  were  taken  into  account. 


Using  the  conservation  of  the  flux  of  wave  energy  in  the  direction  of  waves  propagation 
and  some  other  general  arguments,  the  changes  in  U  along  the  slope  may  be  estimated 
theoretically  (see  example  in  Fig.  20  and  details  in  Voropayev  et  al.,  2003)  and  local  values  of  \P 
in  (8)  and  L  in  (7)  can  be  accurately  calculated  (Fig.  19b).  It  is  also  shown  in  the  mentioned 
paper  that  when  the  wave  intensity  increases,  x  strongly  decreases,  implying  that  for  more 
energetic  ripples  are  formed  faster  and  reach  a  quasi-steady  state  earlier.  This  explains  why  in 
experiments  with  sandy  slopes  the  ripple  formation  zone  forms  a  front  that  separates  ripples  from 
the  smooth  sandy  bottom  near  the  break  point. 

Figure  20.  A  typical  distribution  of  the 
maximum  water  velocity,  U,  along  the 
slope.  Symbols  -  measured  values,  solid 
line  -  theoretical  estimates.  Experimental 
parameters:  0J=  0.4  Hz,  So  =5  cm.  Waves 
break  in  section  #31. 
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At  large  times  t »  x  the  morphology  of  the  ripple  bed  continues  to  change,  although  the 
mean  ripple  characteristics  Lo  and  ho  (at  a  fixed  position  along  the  slope)  do  not  change 
appreciably  after  the  formation  period  2x;  i.e.  the  bottom  topography  never  reaches  an  exact 
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steady  state.  While  the  frequency  and  amplitude  of  paddle  oscillations  remain  constant,  the 
ripples  generated  on  the  bed  demonstrate  large  time  instabilities. 

The  first  type  of  instability  is  related  to  the  periodic  merger  and  split  of  ripples  (Fig.  21). 
This  effect  was  observed  and  explained  for  the  case  of  a  horizontal  sandy  bed  under  spatially 
homogeneous  oscillatory  flow  in  Voropayev  et  al.  (1999). 


◄ - ► 

Figure  21.  Plan  view 
photographs  showing  a  regular  (a) 
system  of  ripples  in  quasi-steady 
state  (a)  and  the  large-time 
instability  of  the  bottom 
topography  and  formation  of  a 
new  ripple  (b).  With  time  the 
newly  formed  ripple  spreads 
along  the  bed  across  the  flow  O5) 
and  initiates  small  mean  drift  of 
the  surrounding  ripples.  Arrows 
show  the  flow  directions. 


This  process  (Fig.  21)  resembles  the  appearance  of  periodic  dislocations  in  an  initially 
regular  ripple  pattern.  The  dislocation  induces  migration  of  ripples  in  the  proximity,  and  with 
time  the  new  or  merged  ripples  spread  across  the  flow.  Because  the  typical  spacing  between 
ripples  decreases/increases  during  ripple  merging/splitting,  the  whole  system  of  ripples  migrates 
with  some  typical  "stochastic"  drift  velocity  U*d.  With  time,  a  new  quasi-steady  state  with 
approximately  the  same  typical  spacing,  Lo,  between  ripples  and  height,  ho,  is  established,  hi 
spatially  homogeneous  oscillatory  flow,  the  migration  velocity  has  no  preferred  direction  and  the 
mean  drift  velocity  as  averaged  over  large  times  tends  to  be  zero.  The  typical  absolute  value  of 
this  velocity  for  ripples  over  a  horizontal  bottom  may  be  estimated  as  (Voropayev  et  al.,  1999) 

Ud*  =  UvF1/2(2.2  -  (US^VSC.  (9) 

The  observations  collected  in  the  model  surf  zone  show  that,  in  addition  to  the  above- 
mentioned  instability,  another  type  of  instability  is  possible  under  progressive  waves  on  sloping 
beaches.  This  instability  causes  a  steady  onshore  ripple  drift  with  unidirectional  drift  velocity, 
Ud,  along  the  slope.  Using  the  model  proposed  by  Kennedy  (1963)  and  Engelung  (1970)  (see 
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also  Nielsen,  1992,  p.  133)  for  sediment  transport  resulting  from  unidirectional  drift  of  sinusoidal 
bedform  with  constant  shape,  the  value  of  Ud  was  estimated  as  (Voropayev  et  al.,  2003) 

Ud  =  md2(xf'/10)3/2/nho  (10) 

(n  is  the  volume  of  solid  fraction  in  a  unit  volume  of  the  bed,  for  sand,  n  *  0.7).  The  data  on  Ud 
were  collected  for  different  along  slope  positions  (section  #),  and  typical  results  for  one 
experiment  are  shown  in  Fig.  22.  The  estimates  based  on  (10)  (where  local  values  of  £  and  U 
were  used  to  calculate  local  'F)  are  shown  by  solid  lines. 


Figure  22.  Unidirectional  ripple  drift 
velocity,  U&  as  a  function  of  the  along 
slope  position.  Symbols  show  measured 
values  for  an  experiment  with  co  -  0.4  Hz, 
e0  =  10  cm,  solid  line  shows  the 
estimation  based  on  (10). 


To  gain  a  better  understanding  of  sediment  transport  and  ripple  drift,  the  flow  around 
ripples  was  investigated.  Recently,  Ridler  &  Sleath  (2000)  reported  observations  on  the  Eulerian 
time-mean  drift  induced  by  waves  over  solid  artificial  ripples  placed  at  the  bottom  of  a  water 
layer  of  constant  depth.  In  their  work,  velocity  measurements  were  taken  at  different  distances 
from  a  crest  and  a  trough  and  then  phase  averaged.  A  mean  negative  (offshore)  water  drift  was 
observed  near  the  bed  over  ripple  troughs.  Just  above  the  ripple  crests  in  the  boundary  layer, 
however,  a  positive  (onshore)  mean  water  drift  was  documented.  No  similar  data  exists  for 
naturally  formed  (moving)  sand  ripples.  Owing  to  the  drift  and  spatial  variability  of  sand  ripples, 
however  these  measurements  are  more  intricate  and  some  results  are  described  below. 

The  ADV  probe  was  placed  in  the  middle  of  the  tank  in  section  #27  (Fig.  23).  In  addition 
to  the  velocity,  the  ADV  probe  can  measure  the  distance  to  the  bottom  with  an  accuracy  better 
than  0.1  cm.  Using  this  latter  capability,  measurements  of  the  ripple  profile  were  taken  every  10 
s  while  ripples  are  drifting  onshore  under  the  stationary  probe  (Fig.  23). 


25 


Figure  23.  Side  view  photograph  ( digitized 
image  from  video)  showing  the  ADV  probe  above 
the  drifting  (from  right  to  left)  sand  ripples. 

Experimental  parameters:  (Q  =  0.4  Hz,  £o  =  10 
cm,  section  #27.  Frame  (b)  was  taken  2.5  min 
after  the  frame  (a).  The  ADV  probe  is  barely 
visible  and  is  shown  by  an  arrow. 

Using  these  data  and  measuring  the  ripple  length,  the  averaged  ripple  drift  velocity  over 
the  measurement  period  was  calculated  as  Ud  =  0.048  cms'1.  After  collecting  data  on  the  ripple 
profile  and  drift,  the  ADV  probe  was  rapidly  traversed  horizontally  so  that  its  measuring  volume 
is  at  a  location  just  above  the  ripple  crest.  Then  the  velocity  measurements  were  taken  while 
ripples  drifted  onshore.  Figure  24  gives  typical  horizontal  velocity  traces.  The  data  shown  in 
Fig.  24a  were  taken  approximately  5  cm  above  the  ripple  trough  at  the  level  of  the  ripple  crest 
and  the  data  shown  in  Fig.  24b  were  taken  in  the  boundary  layer  approximately  just  above  the 
ripple  crest.  The  velocity  magnitude  at  5  cm  above  the  trough  is  approximately  two  times  larger 
than  in  the  boundary  layer  above  the  crest.  The  averaged  velocity  over  a  wave  period  above  the 
trough  (Fig.  24a)  is  negative  (offshore)  while  it  is  strongly  positive  just  above  the  ripple  crest 
(onshore)  (Fig.  24b),  pointing  to  the  possibility  of  significant  onshore  sediment  transport  and 
hence  a  corresponding  ripple  drift. 


Figure  24.  Typical  examples  of 
horizontal  along-slope  velocity  traces  as 
obtained  in  the  boundary  layer  above 
the  ripple  crest  (a)  and  approximately  5 
cm  above  the  ripple  trough  at  the  level 
of  ripple  crest  (b).  The  velocity 
magnitude  above  the  trough  (b)  is 
approximately  twice  larger  than  in  the 
boundary  layer  above  the  crest  (a).  As 
can  be  seen,  averaged  ( over  the  wave 
period)  velocity  above  the  trough  (b)  is 
negative,  while  it  is  strongly  positive  just 
above  the  ripple  crest  (a).  Time  is  given 
in  accordance  with  Fig.  25  below. 
Experimental  conditions  are  the  same  as 
in  Fig.  23. 
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The  data  in  Fig.  24  present  short  pieces  of  a  rather  long  time  record  collected  while  the 
ripple  drifted  approximately  a  distance  Lo  under  the  ADV  probe.  These  data  were  time-averaged 
over  five  wave  periods  and  the  results  are  shown  in  Fig.  25.  Note  the  positive  mean  water  drift 
just  above  the  ripple  crest.  This  result  is  consistent  with  the  measurements  over  solid  ripples  by 
Ridler  &  Sleath  (2000)  and  justifies  the  use  of  (10)  for  estimating  ripple  drift  velocity. 


Figure  25.  Time  record  of  the  along- 
slope  velocity  component,  u,  averaged 
over  five  wave  periods  and  taken  at  a 
level  just  above  the  ripple  crest  when  the 
ripple  slowly  drifted  one  ripple  length  L0 
under  the  ADV  probe.  Experimental 
parameters  are  same  as  in  Fig.  23 


hi  order  to  better  understand  cobble-sandbar  interactions  (which  are  discussed  in  the  next 
section),  experiments  were  conducted  on  the  formation  and  evolution  of  sandbars  on  an  initially 
flat  sandy  slope,  hi  the  range  of  parameters  used  in  experiments  at  least  two  different  types  of 
sandbars  were  observed,  namely  onshore  moving  sandbars  and  offshore  moving  sandbars.  The 
type  of  breaker  on  a  sloping  beach  depends  on  the  beach  slope  and  wave  parameters.  Battjes 
(1974)  characterizes  the  waves  by  the  surf  similarity  parameter 

£)  =  tana/yfr\*JX, 

where  a  is  the  slope  angle,  r\*  wave  amplitude  near  the  break  and  A,  wavelength  in  deep  water. 
Accordingly,  a  spilling  breaker  occurs  when  £  <  0.4,  a  plunging  breaker  occurs  when  0.4  <  £  < 
2.0  and  a  surging  breaker  occurs  at  larger  Depending  on  the  range  of  parameters  used  in  the 
experiments,  spilling  or  plunging  breakers  can  be  generated,  leading  to  different  sandbar  types; 
an  example  of  a  sandbar  under  a  spilling  breaker  is  shown  in  Fig.  26a,  which  has  a  bump-like 
shape.  In  this  case,  even  after  breaking,  waves  are  still  energetic  enough  to  form  a  large  vortex 
onshore  of  the  accumulated  sand,  causing  offshore  sediment  motion  and  slow  offshore  drift  of 
the  bar.  Note  that  this  opposes  the  sediment  transported  by  the  onshore  migrating  ripples 
offshore  of  the  bar.  In  the  case  of  a  plunging  breaker  (Fig.  26b),  sandbars  of  trapezoidal  shape 
are  formed.  They  have  flat  tops,  small  height  and  larger  widths.  Plunging  breakers  form  closer 
to  the  shoreline  and  significant  amount  of  energy  is  dissipated  by  a  large  vortex.  The  flow  is  not 
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energetic  enough  to  form  another  large  vortex  onshore  side  of  the  accumulated  sand.  Therefore, 
onshore  mean  sediment  transport  is  observed  and,  as  result,  bar  slowly  expands  in  the  onshore 
direction. 


Figure  26.  Two  typical  examples  of 
sandbars  formed  under  spilling  (a)  or 
plunging  (b)  breakers.  Experimental 
parameters :  (O-  0.75  Hz,  Co  =  7.5  cm 
(a),  0.6,  5  (b).  Time  is  approximately  t 
=  140  min  (a),  75  (b)  from  the 
beginning  of  the  experiment. 


5.  Cobbles  on  sandy  bed 

Upon  completion  of  detailed  studies  on  the  sand  bottom  morphodynamics,  model 
cobbles/mines  were  placed  at  the  slope  and  their  behavior,  including  scouring  and  burial,  was 
investigated.  The  experiments  were  initiated  with  a  smooth  sandy  bottom.  When  the  system 
reached  its  quasi-steady  state,  with  established  ripples,  cobbles  were  dropped  on  the  bed  and 
their  behavior  was  monitored.  Since  the  flow  velocity  around  the  cobble  increases  relative  to  the 
background  flow,  scouring  of  the  bed  can  be  initiated.  The  results  of  experiments  show  (Fig.  27) 
that  the  maximum  subsidence,  S,  of  the  cobble  is  not  very  significant  and  doesn't  depend  on  the 
cobble  density  and  height,  but  depends  mostly  on  the  cobble  diameter,  D,  and  may  be  estimated 
as  S  =  0.15D  (Voropayev  et  al.,  1999).  At  the  same  time,  sediment  from  the  nearest  ripples  is 
transported  towards  the  cobble,  thus  significantly  changing  the  bottom  topography  close  to  the 
cobble  and  burying  the  cobble  under  migrating  ripples  (Fig.  28).  Thus,  the  scouring  appeared  to 
be  of  lesser  importance  as  compared  to  the  migration  of  ripples. 

This  phenomenon  (ripple  migration)  leads  to  periodic  burial  of  cobbles  by  migrating 
ripples.  Typical  period  of  such  burial  can  be  estimated  as  Tb  =  Lo/Ud,  which  is  the  time  for  a 
system  of  ripples  to  drift  with  a  velocity  Ud  a  distance  equal  to  one  ripple  length  Lo  (see 
schematic  in  Fig.  29).  Estimates  show  that,  for  example,  for  conditions  used  in  the  experiment  in 
Fig,  27,  Tb  =  6.5  min,  which  is  in  agreement  with  the  observations.  For  the  typical  oceanic 
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conditions  (wave  amplitude  1  m,  period  10  s,  depth  10  m,  sand  size  0.1  cm)  estimates  also  show 
(Voropayev  et  al.,  1999)  that  a  mine  of  diameter  20  cm  and  height  10  cm  may  be  completely 
buried  after  a  typical  time  on  the  order  of  4  hours,  which  is  in  agreement  with  the  reported  field 
data  on  ripple  migration  (Traycovski  et  al.,  1999). 


Figure  27.  Front  view  (a,  b) 
and  plan  view  (c)  showing  the 
scour  pattern  evolution  near  the 
cobble  in  steady  oscillatory 
flow.  Three-dimensional  small- 
scale  rolling  grain  ripples  are 
also  visible  away  from  the 
cobble  in  (b).  Arrows  show  the 
directions  of  the  flow. 
Experimental  parameters:  U0  = 
15  cms'1,  d  =  0.04  cm,  D  =  10 
cm,  H  =  6  cm.  pc  =  2.0  gmcm3, 
time  t  =  5  min  (a),  55  min  (b,c) 
after  the  beginning  of  the 
experiment.  The  ADV  probe  is 
also  shown  on  these 
photographs. 


(a) 


Figure.  28.  Periodic  burial  of 
relatively  small,  D  =  8  cm,  cobble 
that  was  trapped  between  two 
adjacent  migrating  ripples.  Time 
interval  between  frames  (a)  and 
(d)  is  6  min. 


(b) 


Figure  29.  Schematic  showing  the  time 
behavior  of  a  cobble/mine,  which  was  placed 
on  a  sandy  bed  in  oscillatory  flow:  (a)  local 
scour  and  subsidence  of  cobble;  (b-d) 
periodical  burial  of  cobble  under  migrating 
sand  ripples. 


Another  topic  that  continues  to  being  addressed  in  our  research  is  sand  bar  interaction 
with  large  solid  objects.  As  mentioned,  at  large  times  a  sand  bar  of  size  much  larger  than  the 
cobble  is  established  under  the  breaking  waves  (see  Fig.  26).  The  water  depth  decreases  rapidly 
along  the  slope  of  the  bar,  and  hence  the  flow  becomes  more  energetic  and  provides  enough 
momentum  for  the  cobbles  to  climb  along  the  sloping  bar  and  reach  its  top. 

If  the  incoming  wave  characteristics  remain  fixed,  a  quasi-equilibrium  state  can  be 
established  during  which  the  cobbles  stay  near  the  top  of  the  bar  for  a  period  of  time  (in 
experiments  0.5  hour  or  more);  they,  however,  oscillate  back  and  forth  under  the  wave  action 
(Fig.  30a).  This  observation  can  be  explained  as  follows.  The  onshore  water  motion  above  the 
bar  is  very  energetic  and  forces  the  cobble  onshore.  In  addition,  the  incoming  waves  break  here 
and  a  large  vortex  is  formed  under  the  wave  crest,  resulting  in  an  intense  turbulent  bore  that 
propagates  onshore  (swash  flow).  Return  flow  in  the  swash  is  less  energetic,  yet  sufficiently 
strong  and  propagates  near  the  bottom  under  the  bore.  The  returning  swash  flow  and  the 
offshore  motion  (generated  under  the  wave  trough)  cause  the  cobble  to  move  offshore,  following 
its  onshore  motion. 

The  oscillations  of  the  cobbles  at  the  top  of  the  bar  occur  by  these  two  competing 
mechanisms.  This  position,  however,  is  unstable.  Due  to  intense  turbulent  motion,  the  cobble 
may  tip  and  slide  down  along  the  bar.  When  the  cobble  slides  to  the  offshore  side  of  the  bar,  it  is 
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usually  pushed  back  to  the  top  of  the  bar  and  start  oscillating  again.  If  it  slides  along  the  onshore 
side  of  the  bar,  the  cobble  can  either  be  pushed  back  to  the  top  of  the  bar  or  get  entrapped  there. 


Figure  30.  Large  and  relatively  light  cobble,  which  was  initially  above  the  sand  bar 
when  the  incoming  wave  characteristics  were  changed.  The  cobble  oscillates  steadily 
above  the  bar  (a),  then  it  slides  to  the  trough  at  the  onshore  side  of  the  bar  (b)  and 
finally  moves  steadily  onshore  (c),  reaching  with  time  the  shoreline  ( not  shown). 
Cobble  2  was  used  in  this  run  (see  Table  1). 


Changing  of  incoming  wave  characteristics  alters  the  position  of  the  break  point  and 
hence  the  behavior  of  cobbles.  The  associated  processes  are  intricate,  and  here  we  present  an 
example  to  illustrate  this  complexity. 


Figure  31.  The  burial  of  a  large  and 
relatively  heavy  cobble,  which  initially 
was  above  the  sand  bar  when  the 
incoming  wave  characteristics  were 
changed.  The  cobble  (a)  slides  to  the 
trough  at  the  onshore  side  of  the  bar 
and  is  slowly  buried  (b)  under  the 
drifting  bar.  Cobble  3  was  used  in  this 
run. 


When  the  wave  characteristics  are  changed  to  a  lower  frequency  or  smaller  amplitude,  the 
waves  break  closer  to  the  coastline  and  the  bar  drifts  onshore.  Thence,  the  back  flow  above  the 
bar  is  significantly  reduced  and  the  cobble  slides  down  from  the  bar  crest  along  the  onshore  side 


31 


(Fig.  30b).  The  general  behavior  of  cobble  depends  on  the  cobble  density.  If  the  cobble  is 
relatively  light,  it  moves  slowly  onshore,  reaches  the  swash  region  and  (usually)  continues  to 
move  (Fig.  30c)  under  the  turbulent  bore  until  it  arrives  at  the  shoreline.  If  the  cobble  is 
relatively  heavy,  it  tends  to  be  trapped  in  the  onshore  face  of  the  bar.  Because  of  the  slow 
onshore  migration  of  the  bar,  the  cobble  can  be  completely  buried  under  the  bar,  as  shown  in  Fig. 
31,  for  a  long  period  of  time.  This  burial  process  is  rather  slow  and  can  take  several  hours  in 
laboratory. 

6.  Conclusions 

The  main  finding  of  this  research  may  be  summarized  briefly  as  follows. 

A.  hi  the  first  phase  of  the  research,  the  motion  of  cobbles/mines  in  a  swash  zone  on  an 
impermeable  slope  was  modeled  experimentally  and  analyzed  theoretically,  hi  oceanic 
situations,  the  beach  is,  of  course,  permeable.  Owing  to  the  fact  that  very  little  research  has  been 
done  on  this  problem,  it  was  decided  to  first  investigate  the  more  simplified  case  of  an 
impermeable  bottom.  This  allowed  considering  a  physical  system  with  a  reduced  number  of 
external  parameters,  yet  one  which  captured  the  main  features  of  oceanic  situations.  Even  with 
this  simplified  problem,  at  least  ten  external  parameters  (H0,D,hc,pc,p,v,g,a,K  and  t)  were 

shown  to  be  important.  In  the  experiments,  a  dam-break  flow  was  used  to  simulate  an  impulsive 
swash.  It  is  shown  that,  when  the  threshold  for  the  initiation  of  motion  is  exceeded,  the  cobble 
starts  to  move  with  relatively  large  velocity  in  the  on-shore  direction.  The  cobble  displacements 
were  measured  as  a  function  of  time  and  all  other  external  parameters.  Based  on  the  results  of 
experiments  a  physical  understanding  was  obtained  and  a  theoretical  model  was  proposed  to 
calculate  cobble  displacements.  Satisfactory  agreement  between  the  laboratory  observations  and 
the  model  predictions  was  obtained. 

B.  The  bottom  motion  of  spherical  and  disk-shaped  cobbles  in  a  wave-induced  oscillatory 
flow  was  studied  experimentally  and  analyzed  theoretically.  In  the  experiments,  large  standing 
waves  were  used  to  model  an  oscillatory  background  flow  that  occurs  near  the  surf  zone.  It  is 
shown  that  the  cobble  motions  in  such  a  flow  rapidly  become  periodic  with  reduced  amplitude 
compared  to  the  background  flow;  a  phase  shift  between  the  cobble  and  background  flow 
displacements  is  also  evident.  The  cobble  displacements,  as  a  function  of  time  and  other  (at  least 
ten)  external  parameters,  were  measured.  Based  on  the  results  of  the  experiments,  a  model  was 
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developed  and  the  cobble  displacements  were  calculated  and  compared  with  measured  values. 
The  model  also  allowed  calculating  the  rather  complex  behavior  of  the  angular  velocity  of  the 
spherical  cobbles.  The  purpose  was  not  only  to  make  measurements  relevant  to  the  dynamics  of 
cobble  in  oscillatory  flow,  but  also  to  understand  the  basic  physics  and  verify  a  theoretical 
model.  To  reduce  the  number  of  governing  parameters,  a  simplified  geometry  with  a  solid 
bottom,  was  chosen  and  a  model  (which  still  includes  at  least  ten  parameters:  D,  pc,  he,  K,  g, 

p,v,(D,e,  t)  was  developed. 

C.  The  flow  in  the  entire  surf  zone  was  modeled  in  a  large  wave  tank  and  the  motion  of 
cobbles/mines  along  an  impermeable  slope  was  studied  under  time  and  spatial  dependent  flow 
conditions.  Onshore  and  offshore  motions  of  cobbles,  as  well  as  steady  oscillations  with  zero 
mean  displacement,  were  observed  for  different  conditions.  To  explain  the  results  of 
observations,  the  models  developed  in  A  and  B  were  generalized  for  these  more  complicated 
conditions.  Proper  parameterizations  were  used  for  the  pressure  accelerating  term,  drag,  lift  and 
other  nonlinear  forces.  Comparison  of  measured  and  calculated  values  for  the  cobble 
displacement  shows  that  the  model  correctly  predicts  the  cobble  behavior  in  model  surf.  One  of 
the  key  variables  in  the  model,  which  is  known  with  the  least  accuracy,  is  the  virtual  mass 
coefficient  of  a  disc  moving  with  variable  velocity  along  a  solid  boundary.  A  fixed  value  for  this 
coefficient,  as  well  as  for  the  drag  and  lift  coefficients  was  used  in  the  calculations.  It  is 
important  that  in  calculations  the  spatial  variation  of  the  background  velocity  along  the  slope 
must  be  taken  into  account.  Proper  parameterization  for  this  effect  was  proposed  and  used. 

D.  The  behavior  of  cobbles/mines  placed  on  a  horizontal  sandy  bed  in  an  oscillatory  flow 
was  studied  experimentally.  It  is  shown  that  when  the  background  velocity  is  below  its  critical 
value,  small-scale  rolling-grain  ripples  develop  on  an  initially  flat  bottom  and  noticeable  scour 
around  the  cobble  occurs.  At  large  times  the  scour  pattern  reaches  a  quasi-steady  state  and  the 
estimates  for  the  resulting  maximum  subsidence  of  the  cobble  are  derived.  When  the 
background  velocity  exceeds  its  critical  value,  vortex  ripples  start  to  form.  The  typical  time  of 
ripple  formation  depends  on  the  wave  frequency  and  the  mobility  parameter,  and  can  be 
estimated  using  proposed  parameterizations.  At  larger  times  the  system  of  ripples  is  not  stable 
and  they  migrate  with  some  drift  velocity.  As  a  result  a  heavy  cobble  may  be  buried 
periodically.  The  data  on  the  ripple  drift  velocity  were  collected  and  parameterized.  Basing  on 
this  parameterization  the  estimates  for  the  time  of  periodic  cobble  burial  were  derived. 
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E.  The  evolution  of  an  initially  flat  sandy  slope  and  the  behavior  of  disk-shaped  cobbles 
in  model  surf  zone  were  studied  experimentally.  The  experiments  were  conducted  in  a  large 
wave  tank  with  a  sandy  slope.  Upon  initiation  of  wave  forcing,  the  initially  flat  topography 
changes  and  after  a  transitional  time,  the  bedform  reaches  a  quasi-steady  state  with  a  system  of 
sand  ripples  and  a  large  bar  near  the  break.  Although  the  incoming  wave  characteristics  are 
fixed,  the  bottom  topography  never  reaches  a  full  steady  state  but  evolves  slowly  (e.g.  ripple  drift 
and  sandbar  transformation  at  large  times).  The  data  on  the  ripple  formation,  growth,  drift,  and 
quasi-steady  ripple  characteristics  were  collected  and  compared  with  proposed  theoretical 
predictions.  When  the  incoming  wave  characteristics  were  changed,  the  bottom  topography 
responded  by  establishing  a  new  quasi-steady  state  after  a  certain  adjustment  time  period.  The 
above  studies  were  extended  to  include  cobbles/mines  placed  along  the  sandy  slope.  Four 
different  scenarios  were  identified  and  explained:  (i)  steady  oscillations  of  cobbles  with  zero 
mean  displacement  and  small  scour;  (ii)  mean  onshore  motion  of  relatively  light  cobbles,  (iii) 
periodic  burial  of  relatively  heavy  cobbles  of  which  the  dimensions  are  comparable  to  that  of 
sand  ripples,  (iv)  burial  of  relatively  large  cobbles  under  the  bar  when  the  break  point  was 
changed  by  changing  incoming  wave  characteristics.  Physical  explanations  were  provided  to 
explain  the  observations.  Data  on  the  flow  characteristics  in  the  vicinity  of  the  movable  bed 
were  also  obtained.  The  presence  of  a  mean  onshore  flow  in  the  boundary  layer  above  the  ripple 
crest  was  observed,  corroborating  recent  observations  of  Ridler  &  Sleath  (2000)  carried  out  in 
the  wave  boundary  layer  above  solid  sinusoidal  bedforms. 
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